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OPERATOR ERROR ESTIMATES FOR 
HOMOGENIZATION OF THE ELLIPTIC DIRICHLET 
PROBLEM IN A BOUNDED DOMAIN 



O M - A - PAKHNIN AND T. A. SUSLINA 

(N ' 

^ ■ Abstract. Let C R d be a bounded domain of class C 2 . In the 

Hilbert space 1/2(0; C n ), we consider a matrix elliptic second order dif- 
ferential operator Ad,e with the Dirichlet boundary condition. Here 
e > is the small parameter. The coefficients of the operator are peri- 
odic and depend on x/e. We find approximation of the operator Ap E 
in the norm of operators acting from 1/2(0; C") to the Sobolev space 
Qh ' H 1 (0;C") with an error term O(^fe). This approximation is given by 

, the sum of the operator (A^)^ 1 and the first order corrector, where A%> 

T ■ is the effective operator with constant coefficients and with the Dirichlet 

l~i ' boundary condition. 



Introduction 

> ; 

The paper concerns homogenization theory of periodic differential oper- 
ators (DO's). A broad literature is devoted to homogenization problems in 
■ the small period limit. First of all, we mention the books [BeLP], [BaPa], 

[ZhKO]. 

' 0.1. Operator-theoretic approach to homogenization problems. 

CN . In a series of papers [BSul-5] by M. Sh. Birman and T. A. Suslina a new 

operator-theoretic (spectral) approach to homogenization problems was sug- 
gested and developed. By this approach, the so-called operator error esti- 
mates in homogenization problems for elliptic DO's were obtained. Matrix 
^ ' elliptic DO's acting in L2(M d ; C n ) and admitting a factorization of the form 

A £ = 6(D)*(7(x/e)6(D), e > 0, were studied. Here g(x) is an (m x m)- 
matrix-valued function; it is assumed to be bounded, uniformly positive 
definite and periodic with respect to some lattice V. By we denote the 
elementary cell of the lattice V. It is assumed that m > n and 6(D) is an 
(m x n)-matrix homogeneous first order DO such that rank6(£) = n for 
7^ £ € M. d . The simplest example of such operator is the scalar elliptic 
operator A £ = — divg(x/e)V. The operator of elasticity theory also can be 
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written in the required form. These and other examples are considered in 
[BSu2] in detail. 

In [BSul-5], the equation A £ u £ + u £ = F, where F <G L 2 (M d ;C n ), was 
considered. The behavior of the solution u £ for small e was studied. The 
solution u £ converges in L2(M d ; C") to the solution uo of the "homogenized" 
equation -4°u + u = F, as e ->■ 0. Here A = 6(D)* g°b(D) is the effective 
operator with the constant effective matrix g°. In [BSul,2], it was proved 
that 

ll u e _ u o|lL 2 (IR d ) ^ C £ ll F llL 2 (IR d )- 

In operator terms it means that the resolvent (A £ + J) -1 converges in the 
operator norm in L2(M d ;C n ) to the resolvent of the effective operator, as 
e — s- 0, and 

||(^ + /)- 1 -(^° + /)- 1 || L2 (r< I )_ L2( r« J ) <Ce. (0.1) 

In [BSu3,4], more accurate approximation of the resolvent (A £ + I) -1 in 
the operator norm in L2(M. d ; C n ) with an error term 0(e 2 ) was obtained. 

In [BSu5], approximation of the resolvent (A £ + I)^ 1 in the norm of op- 
erators acting from L,2(M. d ; C n ) to the Sobolev space H 1 (R d ; C n ) was found: 

\\(A £ + I)" 1 - (A° + I)' 1 - eK(e)\\ L2{Rd) ^ HHRd) < Ce; (0.2) 

this corresponds to approximation of u £ in the "energy" norm. Here K{e) 
is a corrector. It contains rapidly oscillating factors and so depends on e. 

Estimates (0.1), (0.2) are called the operator error estimates. They are 
order-sharp; the constants in estimates are controlled explicitly in terms of 
the problem data. The method of [BSul-5] is based on the scaling transfor- 
mation and the Floquet-Bloch theory. The operator A = b(D)* g(x)b(D) is 
decomposed in the direct integral of the operators »4(k) acting in L2(0;C n ) 
and depending on the parameter k (the quasimomentum) . The operator 
family *4(k) has discrete spectrum and depends on k analytically. It is 
studied by methods of the analytic perturbation theory. It turns out that 
only the spectral characteristics of the operator A near the bottom of its 
spectrum are important for constructing the effective operator and obtain- 
ing error estimates. This shows that homogenization can be studied as a 
spectral threshold effect. 

0.2. A different approach to operator error estimates in homogeniza- 
tion problems was suggested by V. V. Zhikov. In [Zhl, Zh2, ZhPas, Pas], 
the scalar elliptic operator — divg(x/e)V (where g(x) is a matrix with real 
entries) and the system of elasticity theory were studied. Estimates of the 
form (0.1), (0.2) for the corresponding problems in M d were obtained. The 
method was based on analysis of the first order approximation to the so- 
lution and introducing of an additional parameter (the shift by the vector 
cj € Q). Besides the problems in homogenization problems in a bounded 
domain O C M d with the Dirichlet or Neumann boundary condition were 
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studied. Approximation of the solution in H l {0) was deduced from the cor- 
responding result in IR d . Due to the "boundary layer" influence, estimates 
in a bounded domain become worse and the error term is 0(e 1 ^ 2 ). The es- 
timate ||u £ — Uo||l 2 (o) < Ce 1 / 2 ||F|| i2 ( C )) follows from approximation of the 
solution in H l (0) by roughening. 

Similar results for the operator — div<7(x/e)V in a bounded domain with 
the Dirichlet or Neumann boundary condition were obtained in the papers 
[Grl, Gr2] by G. Griso by the "unfolding" method. 

0.3. Main results. In the present paper, we study matrix DO's Ad,s in 
a bounded domain O C M d of class C 2 . The operator Ad,e is defined by the 
differential expression 6(D)*g(x/e)6(D) with the Dirichlet condition on dO. 
The effective operator A® D is given by the expression 6(D)*<7°6(D) with the 
Dirichlet boundary condition. The behavior for small e of the solution u e of 
the equation Ad,e u £ = F, where F € £2(0; C n ), is studied. Estimates for 
the i? 1 -norm of the difference of the solution u £ and its first order approx- 
imation are obtained. Roughening this result, we estimate ||u £ — uo||£ 2 (e>). 
Here Uo is the solution of the equation A^uq = F. 

Main results of the paper are Theorems 6.1 and 7.1. In operator terms, 
the following estimates are obtained: 

\\A D \ - {Al)- 1 - eK D {e)\\ L2i0) ^ m{0) < Ce 1 ' 2 , (0.3) 

\\A- D \ - (A%)- l \\ L2{0) ^ L2{0) < Ce 1 ' 2 . (0.4) 

Here K£>(e) is the corresponding corrector. The form of the corrector de- 
pends on the properties of the periodic solution A(x) of the auxiliary problem 
(1.5). If A is bounded, the corrector has a standard form (Theorem 6.1). 
In the general case the corrector contains an auxiliary smoothing operator 
(Theorem 7.1). Besides approximation of the solution u £ in i7 1 (C;C n ), we 
also obtain approximation of the "flux" p £ := g £ b{Y))\y e in Li{0\ C m ). 

0.4. The method is based on using estimates (0.1), (0.2) for homogeniza- 
tion problem in M rf obtained in [BSu2,5] and on the tricks suggested in [Zh2], 
[ZhPas] that allow one to deduce estimate (0.3) from (0.1), (0.2). Main diffi- 
culties are related to estimating of the " discrepancy" w £ which satisfies the 
equation *4 £ w £ = in O and the boundary condition w £ = eK£,(e)'F on 
dO. Note that we can not use the facts specific for scalar elliptic equations, 
because we study a wide class of matrix elliptic DO's. 

0.5. Error estimates in £2(0). It must be mentioned that estimate (0.4) 
is quite a rough consequence of (0.3). So, the refinement of estimate (0.4) 
is a natural problem. In [ZhPas], for the case of the scalar elliptic operator 
— div(?(x/e)V (where <?(x) is a matrix with real entries) an estimate for 

1 n 1 — 

11*4^, e — (A D )~ \\l 2 ^l 2 °f order e 2d ~ 2 for d > 3 and of order e|log£| for 
d = 2 was obtained. The proof essentially relies on the maximum principle 
which is specific for scalar elliptic equations. 



1 



M. A. PAKHNIN AND T. A. SUSLINA 



Using the results and technique of the present paper, one of the authors 
has obtained a sharp order operator error estimate 



A separate paper [Su] is devoted to the proof of this result. 

0.6. The plan of the paper. The paper contains seven sections. In 
Section 1, the class of operators acting in L2(M. d ;C n ) is introduced, the 
effective operator and the corrector are described, and the needed results 
from [BSu2,5] are formulated. In Section 2, properties of the matrix-valued 
function A are described. In Section 3, we introduce the operator smoothing 
in Steklov's sense and prove one more theorem for homogenization problem 
in M. d . Section 4 contains the statement of the problem in a bounded domain 
and description of the "homogenized" problem. In Section 5, we prove some 
auxiliary statements needed for further investigation. Main results of the 
paper are formulated and proved in Sections 6 and 7. Herewith, in Section 
6 the case where A £ Loo is studied, while in Section 7 the general case is 
considered. 

0.6. Notation. Let S) and $)* be complex separable Hilbert spaces. The 
symbols (•, and || • ||^ stand for the inner product and the norm in fj; the 
symbol || • H^-^* denotes the norm of a linear continuous operator acting 
from S) to f)*. 

The symbols (•, •) and | • | stand for the inner product and the norm in C n ; 
1 = l n is the identity (n x n)-matrix. If a is an (n x n)-matrix, the symbol 
\a\ denotes the norm of the matrix a viewed as a linear operator in C n . We 
use the notation x = (xi, . . . , x^) € M. d , iDj = dj = d/dxj, j = 1, . . . ,d, 
D = — iV = (Di, ... ,-Dd). The L p -classes of C"-valued functions in a 
domain O C M d are denoted by L p (0; C n ), 1 < p < oo. The Sobolev classes 
of <C n -valued functions in a domain O C M. d are denoted by H s (0; C n ). By 
HHP; C n ) we denote the closure of (0; C n ) in H l (Q- C n ). If n = 1, we 
write simply L p (0), H s (0), etc., but sometimes we use such abbreviated 
notation also for spaces of vector- valued or matrix- valued functions. 

51. Homogenization problem for a periodic elliptic operator in 



In this section, we describe the class of matrix elliptic operators under 
consideration and formulate the results on homogenization problem in M. d 
obtained in [BSu2,5]. 

1.1. Lattices in M. d . Let ai,...,a ( jGl' i be the basis in M. d that generates 
the lattice T: 



\\A D ,e - (Ad) \\l 2 (0)^l 2 (0) < Ce. 



L 2 (R d ;C n ) 



d 
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and let be the elementary cell of T: 

d ^ 1 

Q := {x € R d : x = Tj aj, ~^<r 3 < -}■ 

3=1 

We denote = measfi. 

The basis bi, . . . , in M. d dual to ai, . . . , is defined by the relations 
(bj,aj) = 2ir5ij. This basis generates the lattice F dual to F: 

d 

r = {b€R rf : b = ^p i b i , Pi €Z}. 
1=1 

We introduce the central Brillouin zone 

f2 = {keM d : |k|<|k-b|, 0/bef}, 

which is a fundamental domain of F. 

Below, H l {£l) stands for the subspace of all functions in H l (Q) whose F- 
periodic extension to R. d belongs to H^ oc (R d ). If <£>(x) is a T-periodic function 
in M. d , we denote 

99 £ (x) := </?(£ _1 x), e > 0. 

1.2. The class of operators. In L2(R d ; C n ), we consider a second order 
DO A e formally given by the differential expression 

A £ = 6(D)*cf(x)&(D), e>0. (1.1) 

Here <?(x) is a measurable (m x m)-matrix-valued function (in general, with 
complex entries). It is assumed that <?(x) is periodic with respect to the 
lattice T, bounded and uniformly positive definite. Next, 6(D) is a homoge- 
neous (m x n)-matrix first order DO with constant coefficients: 

d 

6(D) = (1-2) 

i=i 

Here bi are constant matrices (in general, with complex entries). The symbol 
H€) = Yli=i £ e * s associated with the operator b(F)). We assume 
that m > n and that rank&(£) = n, V£ ^ 0. This is equivalent to the 
inequalities 

a l n < b{6fb{6) < ai l n , 6 e S d_1 , < q < a x < oo, (1.3) 

with some positive constants «o and cti- 

The precise definition of the operator A e is given in terms of the corre- 
sponding quadratic form 

a £ [u,u]= [ ( 5 £ (x)6(D)u,6(D)u) dx, u G H 1 (R d ; C n ). 



6 M. A. PAKHNIN AND T. A. SUSLINA 

Under the above assumptions this form is closed in L,2(M. d ; C ra ) and non- 
negative. Using the Fourier transformation and (1.3), it is easy to check 
that 

c J |Du| 2 dx < o e [u,u] < ci y |Du| 2 dx, u e H 1 (R d ; C n ) , (1.4) 

where c = a o ||0 _1 |lzi> c x = ailbllz^- 

The simplest example of the operator (1.1) is the scalar elliptic operator 
A e = — divg e (x)V = D*g e (x)D. In this case we have n = 1, m = d, 6(D) = 
D. Obviously, (1.3) is true with ao = ot\ = 1. The operator of elasticity 
theory can be also written in the form (1.1) with n = d, m = did + l)/2. 
These and other examples are considered in [BSu2] in detail. 

1.3. The effective operator. In order to formulate the results, we need 
to introduce the effective operator A . 

Let an (n x m)-matrix-valued function A(x) be the (weak) T-periodic 
solution of the problem 

6(D)*0(x) (6(D)A(x) + l ro ) = 0, / A(x) dx = 0. (1.5) 

Jn 

In other words, for the columns Vj(x), j = 1, . . . , m, of the matrix A(x) the 
following is true: Vj G iJ^jC"), we have 

/ ( 5 (x)(6(D)v j (x) + e J ),6(D)r ? (x))dx = 0, V v € H 1 (tyC), 
Jn 

and J a Vj(x) cix = 0. Here ei, . . . , e m is the standard orthonormal basis in 

The so-called effective matrix g° of size mxmis defined as follows: 

g° = l^- 1 I g(x) (6(D)A(x) + l m ) dx. (1.6) 

It turns out that the matrix g° is positive definite. The effective operator 
A for the operator (1.1) is given by the differential expression 

A = 6(D)* 5 °6(D) 

on the domain H 2 (R d ;C n ). 

1.4. Properties of the effective matrix. The following properties of the 
effective matrix are proved in [BSu2, Chapter 3, Theorem 1.5]. 

Proposition 1.1. The effective matrix g° satisfies the following estimates: 

g<g°<g- (1.7) 

Here 

g=\n\~ 1 f g( x )dx, g=(\n\- 1 f gto-Ux) . 
Jn \ Jn J 

If m = n, then g° coincides with g. 
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In homogenization theory, estimates (1.7) are well known for specific DO's 
as the Voight-Reuss bracketing. Now we distinguish the cases where one of 
the inequalities in (1.7) becomes an identity. The following statements were 
checked in [BSu2, Chapter 3, Propositions 1.6 and 1.7]. 

Proposition 1.2. The identity g° = g is equivalent to the relations 

6(D)*g fc (x) = 0, k = l,...,m, (1.8) 

where gfc(x), k = 1, . . . ,m, are the columns of the matrix g(x). 
Proposition 1.3. The identity g° = g is equivalent to the representations 

l fc (x) = lg + 6(D)w fc , \° k eC m , w k eH\n;C n ), k = l,...,m, (1.9) 

where lfc(x), k = 1, . . . ,m, are the columns of the matrix (/(x) -1 . 

Obviously, (1.7) implies the following estimates for the norms of the ma- 
trices g° and (<7°) _1 : 

l/I^NI^, ICs )" 1 ! < llr'llLco- (i-io) 

1.5. The smoothing operator. We need an auxiliary smoothing operator 
n e acting in L2(JR d ;C m ) and defined by the relation 

(n £ u) (x) = (2ny d / 2 f e^u(t) d£, (1.11) 

Jn/e 

where u(£) is the Fourier- image of u(x). In other words, II £ is the pseudodif- 
ferential operator with the symbol Xo/ £ (£) which is the indicator of the set 

tl/e. Obviously, II £ is the orthogonal projection in each space H s (R d ; C m ), 
s > 0. Besides, D a U £ u = U £ D a u for u G H s (R d ;C m ) and any multiindex 
a such that \a\ < s. 

Proposition 1.4. For any u G H 1 (R d ;C m ) we have 

||II £ u - u|| L2( - R d. C m) < er 1 ||Du|| i2 ( M d), 

where ro is the radius of the ball inscribed in clos CI. 
Proof. For £ G R d \ {Cl/e) we have |£| > roe" 1 . Hence, 

R d \(Q/e) 

<e 2 r 2 [ \t\ 2 \u(t)\ 2 dZ = e 2 r 2 [ |Du(x)| 2 dx. . 

Jm d Jm d 

The following property was proved in [BSu5, Subsection 10.2]. 

Proposition 1.5. Let /(x) be a T-periodic function in M. d such that 
f G L2(H). Let [f 6 ] denote the operator of multiplication by the function 
/(£ _1 x). Then the operator [/ £ ]n £ is continuous in L2(M d ;C m ), and 

ll[/ £ ] n e|lL 2 (IR d ;C m )^L2(R d ;C m ) < |^|~ 1/2 ||/|U 2 (n) • 



irLu — u 



|2 

lL 2 (IR d ;C m ) 
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1.6. Results for homogenization problem in M. d . Consider the follow- 
ing elliptic equation in R d : 

A £ u £ + u £ = F, (1.12) 

where F € L2(R d ; C"). It is known that, as e —> 0, the solution u £ converges 
in L2(M d ;C n ) to the solution of the "homogenized" equation 

^°u + u = F. (1.13) 
The following result was obtained in [BSu2, Chapter 4, Theorem 2.1]. 
Theorem 1.6. Let u £ be the solution of the equation (1.12), and let uq be 
the solution of the equation (1.13). Then 

ll u e ~ u ollL 2 (R d ;0) ^ Cl e ll F llL 2 (K d ;C™)) < £ < 1, 

or, in operator terms, 

\\(Ae + I)' 1 - (A + /)" 1 || i2( Rc i;C „ ) ^ L2(M c i . C n ) < Cie, < e < 1. 

The constant C\ depends only on the norms H^Hl^, Ib" 1 !!^^, the constants 
&o, oil from (1.3), and the parameters of the lattice T. 

In order to find approximation of the solution u £ in i7 1 (M d ;C n ), it is 
necessary to take the fist order corrector into account. We put 

K(e) = [A £ ]U £ b(B)(A° + J)" 1 . (1.14) 

Here [A £ ] is the operator of multiplication by the matrix-valued function 
A(e _1 x), and II £ is the smoothing operator defined by (1.11). The opera- 
tor (1.14) is continuous from L 2 (R d ;C n ) to H 1 ^; C n ). This fact can be 
easily checked by using Proposition 1.5 and relation A G H 1 (Q). Herewith, 
e\\K(e)\\ L2 ^ H i=0(l). 

The " first order approximation" of the solution u £ is given by 

v £ = u + eA £ n e 6(D)u = (A + /) _1 F + eK(e)F. (1.15) 

The following theorem was obtained in [BSu5, Theorem 10.6]. 
Theorem 1.7. Let u £ be the solution of the equation (1.12), and let uo be 
the solution of the equation (1.13). Let v £ be the function defined by (1.15). 
Then 

||Ue - v e\\H 1 («. d ;C n ) ^ C2£||F|| L2 ( M d ;C n) , < £ < 1, (1-16) 

or, in operator terms, 

l^ + I)- 1 -(A^I)- 1 -eK(e)\\ L2(Rd . Cn) ^ m(Rd . Cn) < C 2 e, < e < 1. 

The constant C 2 depends only on m,ao,oci, H^Hl^, Ib -1 !!/^, and the pa- 
rameters of the lattice V. 

Now we distinguish the case where the corrector is equal to zero. Next 
statement follows from Theorem 1.7, Proposition 1.2 and equation (1.5). 
Proposition 1.8. If g° = ~g, i. e., if relations (1.8) are satisfied, then A = 
and K{e) = 0. Then we have 

ll U £ _ U 0||_ffl(Rd;C") < C2£||F|| L2 ( R d. C n), o < e < 1. 
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It turns out that under some assumptions on the solution of the problem 
(1.5) the smoothing operator II e in the corrector (1.14) can be removed 
(replaced by the identity). 

Condition 1.9. Suppose that the T-periodic solution A(x) of the problem 
(1.5) is bounded: A € L^. 

We put 

K°(e) = [A £ ]b(V){A° + 

In [BSu5], it was shown that under Condition 1.9 the operator K°(e) is 
a continuous mapping of L2(M d ;C n ) into H 1 (M. d ;C n ). (It is also easy to 
deduce this fact from Corollary 2.4 proved below.) 

Instead of (1.15), we consider another approximation of the solution u e : 

v e = u + eA e 6(D)u = (A° + /) _1 F + eK°(e)F. (1.17) 

The following result was obtained in [BSu5, Theorem 10.8]. 

Theorem 1.10. Suppose that Condition 1.9 is satisfied. Let u £ be the 
solution of the equation (1.12), and let uo be the solution of the equation 
(1.13). Let v e be the function defined by (1-17). Then we have 

\\ u e — V £ ||/fi(]Rd ; c™) — C3£||F|| i2 ( R d. C n), < £ < 1, 

or, in operator terms, 

||(^ £ + /)- 1 -(^ + /)- 1 -^ (£)|| L2 ( ffid;C «Hifi( K rf;C") < C 3 e, 0<e<l. 

The constant C3 depends only on m,d,ao,a\, H^Hl^, ||5 _1 |Uoo) the param- 
eters of the lattice T, and the norm co- 
in some cases Condition 1.9 is valid automatically. The following state- 
ment was checked in [BSu5, Lemma 8.7]. 

Proposition 1.11. Condition 1.9 is a fortiori valid if at least one of the 
following assumptions is satisfied: 

1°. dimension does not exceed two, i. e. d < 2; 

2°. the operator acts in L2(M d ), d > 1, and has the form A e = D*g e (x)D, 
where g(x) is a matrix with real entries; 

3°. dimension is arbitrary and g° = g, i. e., relations (1.9) are satisfied. 

Note that Condition 1.9 can be also ensured by the assumption that the 
matrix g(x) is sufficiently smooth. 

§2. Properties of the matrix- valued function A 

The following statement is proved by analogy with the proof of Lemma 
8.3 from [BSu5]. 

Lemma 2.1. Let A(x) be the T-periodic solution of the problem (1.5). Then 
for any function u G C^°(U d ) we have 

[ |DA(x)| 2 H 2 dx</3 1 || U ||2 +/3 2 f |A(x)| 2 |Dn| 2 dx. (2.1) 
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The constants j3\ and @2 are defined below in (2.12) and depend only on m, 
d, a , ai, Iblkoo, and Wg^U^- 

Proof. Let Vj(x), j = 1, . . . ,m, be the columns of the matrix A(x). By 
(1.5), for any function rj G H l (R d ; C n ) such that r/(x) = for |x| > r (with 
some r > 0) we have 

/ ( 5 (x)(6(D)v j (x) + e j ),6(D)T ? (x))dx = 0. (2.2) 



(2.3) 



Let u G Cg°(M d ). We put »7(x) = Vi (x)|u(x)| 2 . By (1.2), 

d 

6(D)r7(x) = (6(D) Vj (x)) |^(x)| 2 + £ 6 jVj (x)A|«(x) 



Substituting (2.3) in (2.2), we obtain 

f (<?(x) (6(D) Vj (x) + ej) ,6(D) Vj -(x)> |u| 2 dx 

+ / V(g(x) (6(D)vj(x) + ej),6iVj(x)> + ^x = 0. 



Hence, 



J : = f L 1 / 2 6(D)v j 2 |n| 2 dx = - / (g^e^g^b^Vj) \u\ 2 dx 

JR d I JR d ^ ' ' 

" J Rd ^(5 1/2 &(D)v,,5 1/2 fe/v i ) (£),«« + ^x (2 4) 

-/ y^(gej,biVj) (D[uu + uDiu) dx. 

Denote the summands on the right by J±, J2, J3. The first term J\ can be 
estimated as follows: 

|Jx|< / (W + j ^/2 6(D)V 

jR d V 4 



Next, from (1.3) it follows that 

I/ 2 



M 2 dx < \\g\\Lj\u\\ 2 L2(Wid) + -J. 

(2.5) 

(2.6) 



N<«r, i = l,...,d. 
Taking (2.6) into account, we estimate the second term J 2 : 

|«/ 2 | < 2^ \g 1 ' 2 6(D) Vj | |u| ^ l^feiVj-l dx 



1 

4' 



(2.7) 



< ^J + 4dai|M| Loo / |v/|Du| 2 dx. 
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Finally, the term J3 satisfies the estimate 

d 



(2.8) 



\J 3 \<2^\ge j \ \u\ \J2\btVjWDm\jdK 

< HflUoclKIlL^) + d «ill5lUoc J \vj\ 2 \T>u\ 2 dx. 

Combining (2.4), (2.5), (2.7), and (2.8), we obtain 

\j < 2\\g\\ Loo \\u\\l 2{Rd} + 5dai||<7|| Loo J |v/|Du| 2 dx. (2.9) 

Now, we show how the required estimate can be deduced from (2.9). By 
the Fourier transformation, it follows from the lower inequality (1.3) that 

/ \T>{v 3 u)\ 2 dx<a Q l f |6(D)(v^)| 2 dx. 

JR d JR d 

By (1.2), 

d 

b(p)(v jU ) = (b(T>)vj)u + kvjDiu. 

1=1 

Then, taking (2.6) and the expression for J (see (2.4)) into account, we have 
/ \B(v jU )\ 2 dx < 2a Q l / \b(D)vj \ 2 \u\ 2 dx + 2a Q 1 a 1 d / |vj| 2 |Du| 2 dx 

JR d ' JR d ' JR d 

<2a Q 1 \\g- 1 \\ Loo J + 2ctQ 1 a 1 d / |vj [ 2 |Du| 2 dx. 

jR d 

(2.10) 

Obviously, 

/ \~Dv j \ 2 \u\ 2 dx<2 [ |D( Vi n)| 2 dx + 2 f | Vi | 2 |Du| 2 dx. (2.11) 

JR d JR d ' JR d 

Relations (2.9)-(2.11) imply that 

J |Dv/ |u| 2 dx < lea^ 1 !^- 1 !!^!^!!^!^!! 2 ^^) 

+ 2(l + 2dQ^ 1 ai + 20dao 1 ai||^ 1 || Loo ||(7|| Loo ) [ | Vj | 2 |Du| 2 dx. 

JR d 

Summing up over j, we arrive at estimate (2.1) with 
f3i = 16ma 1 \\g- 1 \\ Loo \\g\\ Loc , 

/3 2 = 2(l + 2da " 1 a 1 + 20da - 1 a 1 ||^ 1 || Loo |b|| Loo ). • 
Corollary 2.2. Under Condition 1.9 for any u G i^ 1 (IR c( ) we have 

[ |DA(x)| 2 |u| 2 dx ^PiWuf^ + ^2||A||L / |D U | 2 dx. (2.13) 

JR d v ' JR d 

Proof. Indeed, the second integral in the right-hand side of (2.1) can be 
estimated by ||A||| f Rd |Du| 2 dx. Then (2.13) is valid for any u £ Cq 
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By continuity inequality (2.13) is extended from the dense set C^°(IR <i ) to 
the whole H x (R d ). • 

Next statement follows from Lemma 2.1 by the scaling transformation. 

Lemma 2.3. Under the assumptions of Lemma 2.1 we have 

[ |(DA) £ (x)| 2 Kx)| 2 dx<^||u||2 , Rd I |A £ (x)| 2 |Du| 2 dx. 

JR d ' ' JR d 

Proof. By the changes y = e _1 x and u(x) = v(y), from (2.1) it follows 
that 

/ |(DA)(^ 1 x)| 2 |n(x)| 2 dx= / |(DA)(y)| 2 Ky)|Vd y 

JR d JR d 

<fa[ \v(y)\ 2 e d dy + /3 2 [ |A(y)| 2 |D w «(y)| V dy 

JR d JR d 

= Pi [ |u(x)| 2 dx + /3 2 e 2 f |A(^ 1 x)| 2 |D x u(x)| 2 dx. • 

JM. d Jm. d 

Corollary 2.4. Under Condition 1.9 for any u € ^(W 1 ) we have 

[ |(DA) e (x)| 2 H 2 dx<A|| U ||2 Rd +^ 2 ||A||i oo e 2 / |Dn| 2 dx. 

JR d v ' JR d 

In conclusion of this section, we give two estimates for the matrix-valued 
function A obtained in [BSu4, (6.28) and Subsection 7.3]: 

l|A|| M „) < l^r / W/ 2 (2r )- 1 « - 1 / 2 || 5 ||^|| 5 - 1 ||^, (2.14) 

l|DA|| L2(Q) < iniVW/V 72 !^!!^^" 1 "^- ( 2 - 15 ) 

§3. Smoothing in Steklov's sense. One more result for 
homogenization problem in r d 

In [Zh2, ZhPas], smoothing in Steklov's sense was used instead of the 
smoothing operator (1.11). It turns out that smoothing in Steklov's sense 
is more convenient for the study of homogenization problem in a bounded 
domain. In this section, we show that for the problem in W d both variants 
are possible, i. e., Theorem 1.7 remains true if in the corrector (1.14) the 
operator LI £ is replaced by the operator smoothing in Steklov's sense. 
3.1. Smoothing in Steklov's sense. In L2(M. d ;C m ), we consider the 
operator S £ defined by 

(5 e u)(x) = / u(x-ez)dz (3.1) 

Jn 

and called the operator smoothing in Steklov's sense. It is easy to check that 
\\Seh 2 (R d )^L 2 (R d ) < 1- Obviously, D a S £ u = S £ D a u for u G H s (R d ;C m ) 
and any multiindex a such that \a\ < s. 

We need some properties of the operator (3.1), cf. [ZhPas, Lemmas 1.1 
and 1.2]. 
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Proposition 3.1. For any u € H 1 (M d ;C m ) we have 

\\ S eU ~ u|| L2 ( K d ;C m) < eri||Du|| L2(R d), (3.2) 
where 2r\ = diamfl 
Proof. By the Cauchy inequality, 



in _ ||2 

p £ u - u|| i2( - Md . cm ) 



(3.3) 



= / dx \Q\ 1 / (u(x — ez) — u(x)) dz 

jR d Jn 

< [ dx [ |u(x-ez) -u(x)| 2 dz. 

Using the Fourier transformation, we obtain 

( |u(x-ez) -u(x)| 2 dx = ( |exp(-ie(z,£))-l| 2 |u(£)| 2 d£ 

<e 2 \z\ 2 [ |£| 2 |u(£)| 2 d£ = e 2 |z| 2 f |Du(x)| 2 dx. 

JR d JR d 

Integrating this inequality over z € fi, we conclude that 

[ dz [ |u(x-ez) -u(x)| 2 dx < e 2 r 2 \n\ [ |Du(x)| 2 dx. 

Jn JR d JR d 

Together with (3.3) this implies (3.2). • 

Proposition 3.2. Let /(x) be a T-periodic function in M. d such that f G 
L2(Q). Then the operator [f £ ]S £ is continuous in L2(M d ;C m ), and 

\\[f £ ]Se\\L 2 (R d ;C m )^L 2 (R d ;C m ) ^ 1^1 1 ^ \\ f \\ L 2 {Si) ■ 

Proof. By the Cauchy inequality and the change of variables, from (3.1) it 
follows that 

/ |,f (x)(S £ u)(x)| 2 dx< \Sl\~ 1 f dxl/^x)) 2 / |u(x-ez)| 2 dz 
jR d JR d Jn 

= l^l- 1 £ dy ^ |/(e-V + z)| 2 |u(y)| 2 dz = l^r 1 !!/!! 2 ^^ ||u|| 2 L2(Rd) . . 
3.2. We put 

K(e) = [A £ }S £ b(D)(A° + I)- 1 . (3.4) 

The operator (3.4) is continuous from L 2 {R d ;C n ) to H 1 (R d ;C n ). In- 
deed, the operator b(D)(A° + J) -1 is a continuous mapping of L2(M. d ;C n ) 
into # 1 (M d ;C m ). Using Proposition 3.2 and relation A G ff^fi), it is 
easy to check that the operator [A 5 ]^ is continuous from H 1 (R d ;C m ) to 
i? 1 (R <i ; C n ). 

Let u £ be the solution of the equation (1.12). Instead of (1.15) we consider 
another fist order approximation of u e : 

v £ = u + eA £ ,S e 6(D)uo = (A + /) _1 F + eK(e)F. (3.5) 

Along with Theorem 1.7, the following result takes place. 
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Theorem 3.3. Let u £ be the solution of the equation (1.12), and let uq be 
the solution of the equation (1.13). Let v £ be the function defined by (3.5). 
Then 

ll U £ _ V e || H l( R d. C n) < C2£||F|| L2 ( R d. C n), < £ < 1, (3.6) 

or, in operator terms, 

IKAe + I)- 1 - (A° + I)' 1 - eK(e)\\ L2{Rd . Cn) ^ HHRd . Cn) <C 2 e, 0<e< 1. 

The constant C2 depends only on m, d, ao, cxi, HpIIloo, ||5 ,_1 |!l oo; and the 
parameters of the lattice T. 

Theorem 3.3 will be deduced from Theorem 1.7. 

Lemma 3.4. For any u G H 2 (M. d ;C n ) we have 

f |(DA) £ | 2 |(n £ -S' e )6(D)u| 2 dx</3 1 f |(n £ -S £ )6(D)u| 2 dx 

d (3.7) 
+ /3 2 e 2 V / |A £ | 2 |(n £ -5 £ )6(D)^-u| 2 dx. 

Proof. By Propositions 1.5 and 3.2 and relation A G H l (VL), all the 
terms in (3.7) are continuous functionals of u in the norm of H 2 (R d ; C n ). 
Since <7 °°(R d ;C n ) is dense in H 2 (R d ; C n ), it suffices to check (3.7) for 
u G C^(R d ;C n ). 

We fix a function ( G C 00 ^) such that < ((t) < 1, ((t) = 1 for 
< t < 1, and C(i) = for t > 2. We put Cr(x) = C(^ _1 |x|), x G R d , 
R>0. Let u G Cg°(R d ;C n ). Then 0?(n £ - S" £ )6(D)u G C °°(M d ;C m ) and, 
by Lemma 2.3, we have 

/ |(DA) £ | 2 |C R (n £ - S £ )b(D)u\ 2 dx < ft / |Cij(n £ - 5 £ )6(D)u| 2 dx 

Take into account that max |<9jCr| < cR^ 1 . Then (3.7) follows from the last 
inequality by the limit procedure as R — > 00, by the Lebesgue Theorem. • 
From Proposition 1.5 and estimate (2.14) it follows that 

ll[^- e ]n £ ||L 2 (]Rd ; c™)^L 2 (Rd;C") ^ N 1 ^ 2 II^-IIl 2 (^) 

s 1/2 fr. \-l —1/2 1| M 1/2 I, _l 11 1/2 ,, 

<m' (2r ) a \\g\\£j\g H/^ =: M. 
Similarly, Proposition 3.2 implies that 

ll[^ £ ]5 , e|lL 2 (IK d ;C m )^L2(IK d ;C™) < (3-9) 

Lemma 3.5. M^e Ziaue 

||eA e (n £ - 5 , £ )6(D)uo|| if i( M d ;C n) < ^||uo||fl2(R««;C»)- ( 3 - 10 ) 



(3.8) 
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The constant C is defined below in (3.16) and depends only on m, d, H^Hl^, 
ll<7 -1 Ikoo ) «0) ot-x-, and the parameters of the lattice F. 
Proof. From (1.3), (3.8), and (3.9) it follows that 

||eA £ (n e - S e )6(D)uo|| L2 (Rd ; c») < 2Mai /2 e||u || // i (Rd;C n ) . (3.11) 
Consider the derivatives 

^-(eA £ (n £ - ,S £ )6(D)uo) = (J^-^ (U £ - S £ )b(B)u 



+ eA £ (U £ - 5 e )6(D)a,-uo, j = l,...,d. 



Then 



^2\\d 3 (eA%U £ - S £ )b(V)uo)\\l 2{Rd) 

3=1 

<2 [ |(DA) e | 2 |(n £ - S £ )b(D)u \ 2 dx (3.12) 

d ,. 

+ 2e 2 V / |A e (n e - S , £ )6(D)a i u f 2 dx. 

The second summand in the right-hand side of (3.12) is estimated by using 
(1.3), (3.8), and (3.9): 

d 

2e 2 V / |A e (n e - < S e )6(D)a j uo| 2 dx<8e 2 M 2 a 1 ||u ||^ 2md . c „ v (3.13) 

The first summand in the right-hand side of (3.12) is estimated with the 
help of Lemma 3.4: 

2/ |(DA) e | 2 |(n £ -5 e )6(D)u | 2 dx<2/3 1 f \(U £ - S £ )b(B)u \ 2 dx 

JR d JR d 
d „ 

+ 2f3 2 e 2 y / [A £ | 2 |(n e - S £ )b(D)djU \ 2 dx. 
^jR d 

(3-14) 

Next, by Propositions 1.4 and 3.1 and relation (1.3), we have 

||(n e - 5 e )6(D)uo||L a(R d) < e(r„ 1 + ri)a} /2 ||uo||fia(R<»;C»)- ( 3 - 15 ) 

The second summand in the right-hand side of (3.14) is estimated with the 
help of (3.13). Finally, combining (3.12)~(3.15), we obtain 

d 

^||^( £ A £ (n £ -5 £ )6(D)uo)||i 2(Rd) 
3=1 

< e 2 (8M 2 (1 + f3 2 ) + 2/3 1 (r - 1 + r x ) 2 ) ajuoll 2 ^.^. 
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Together with (3.11) this implies (3.10) with the constant 

C , = a} /2 (M 2 (8/3 2 + 12)+2/? 1 (r - 1 +r 1 ) 2 ) 1/2 . . (3.16) 

Now it is easy to complete the proof of Theorem 3.3. By (1.3) and 
(1.10), we obtain the following lower estimate for the symbol of the effective 
operator: 

KOVtyO > co|*| 2 ln, c = a ||^ 1 !l L i- (3.17) 

Using the Fourier transformation and (3.17), we estimate the norm of the 
function u = (A + I^F in H 2 (R d ; C n ): 

2 



luoHWc*) = / U + l£| 2 ) 2 (HZ) *9°Kt) + In)" 1 ^) 



< / (1 + |^| 2 ) 2 (co|€| 2 + 1)- 2 |F(^)| 2 d€ < (1 + Co 1 ) 2 ||F|li 2(Rd;Cn 
Combining this with (1.15), (3.5), and (3.10), we obtain 

ll v £ ~ v e||Hi(Rd ; c") ^ 1 1 u o 1 1 (Rrf ; c" ) < (1 + c 1 )C , e||F|| L2 ( R ,i ;C n). (3.18) 
Relations (1.16) and (3.18) imply (3.6) with C 2 = C 2 + (1 + c Q l )C. • 

§4. HOMOGENIZATION OF THE DlRICHLET PROBLEM IN A BOUNDED 
DOMAIN: PRELIMINARIES 

4.1. Statement of the problem. Let O C M. d be a bounded domain of 
class C 2 . In L2(C;C n ), we consider the operator Ad,s formally given by 
the differential expression fe(D)*g £ (x)6(D) with the Dirichlet condition on 
dO. Precisely, Ad, £ is the selfadjoint operator in L2(0;C n ) generated by 
the quadratic form 

a D , £ [u,u] = [ (/ (x)6(D)u,6(D)u) dx, u G H^O;C n ). 
Jo 

This form is closed and positive definite. Indeed, let us extend u by zero to 
R d \ O. Then u G H 1 (R d ; C n ). Applying (1.4), we obtain 

c J |Du| 2 dx < o D , e [u,u] < a J |Du| 2 dx, u € -ffo(C;C n ). (4.1) 
o o 
It remains to note that the functional ||Du||£, 2 (<m determines the norm in 
Hl{0;C n ) equivalent to the standard one. 

Our goal is to find approximation for small e for the operator A~^ e in 
the norm of operators acting from L2(0;C n ) to H l {0]C n ). In terms of 
solutions, we are interested in the behavior of the generalized solution u e G 
H^(0;C n ) of the Dirichlet problem 

6(D)V(x)6(D)u e (x) = F(x), xeO; u £ \ do = 0, (4.2) 
where F G L 2 (0; C n ). Then u e = A^\F. 
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4.2. The energy inequality. Now, we consider the problem (4.2) with the 
right-hand side of class C n ) and prove the energy inequality. Recall 

that i? _1 (C;C n ) is defined as the space dual to -£fg(0;C n ) with respect 
to the L 2 (0;C n )-coupling. If f e H-\0;C n ) and rj G H^(0;C n ), the 
symbol f {f,rj) <ix stands for the value of the functional f on the element 
rj. Herewith, 



(f , rj) dx 

o 



< l|f|!//- 1 (0;C' i )ll r ?ll// 1 (e , ; ( C n )- (4.3) 



Lemma 4.1. Let f € tf -1 (0; C n ), and let z £ £ H^(0; C n ) be the generalized 
solution of the Dirichlet problem 

6(D)V(x)6(D)z e (x) = f(x), xeO; z £ \ do = 0. 

In other words, z £ satisfies the identity 

[ (/(x)&(D)z £ ,6(D)77)dx= / {f,r,)dx, V r) G Hl{0;C n ). (4.4) 
Jo Jo 

Then the following estimate called the "energy inequality" is true: 

llZellif^OjC") < C\\{ ||/r-i(0 ; C")- (4-5) 

Fere 0= (1 + (diam O) 2 )^ 1 ll^ 1 || Loo . 
Proof. By the lower estimate (4.1), we have 

||Dz e ||| 2(0) < Cq 1 (/6(D)z £ ,6(D)z £ ) L2(0) . (4.6) 

Next, from (4.3) and (4.4) with r\ = z £ it follows that 

(/6(D)z £ ,6(D)z £ ) 

By the Friedrichs inequality, 

||z £ ||l 2 (o) < (diamO)||Dz e || £a(0) . (4.8) 
Finally, combining relations (4.6)— (4.8), we obtain 
KHhi(O) < (l + (diamO) 2 )||Dz £ ||| 2(0) 

< (1 + (diamO) 2 )c 1 ||f|| H -i (c , ) ||z £ ||^i (cl) . 

This implies (4.5). • 

Roughening the result of Lemma 4.1, we arrive at the following corollary. 
Corollary 4.2. The operator A~^ £ is continuous from L2(0;C n ) to 
H^(0;C n ), and 

\\-A D ]e\\L2(0;C n )->H 1 (0;C n ) < & ■ 

In what follows, we shall need the next statement which is proved with 
the help of Lemma 4.1. 
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Lemma 4.3. Let if) G iT^OjC"), and let r £ G H 1 (0;C n ) be the generalized 
solution of the problem 

6(D)V(x)6(D)r e (x) = 0, xeO; r £ \ do = if)\ d0 . (4.9) 

JTterc 

IKIIffHOiC") < 7blM|jri(©;C»)> 7o = 1 + Cd 1/2 ai||5|| Loo . (4.10) 
Proof. By (4.9), the function r £ — if) is the solution of the Dirichlet problem 
A £ (r £ - if)) = -A £ if) mO; (r £ - if))\ 90 = 0. (4.11) 
Here the right-hand side in the equation belongs to C n ), and 

(g*b(D)i{),b(p)ip) L2{0) 



(4.13) 



\\A £ il)\\ H -i {0 ) = sup 

0^6Hj(O;C») IIVllfTi(O) (4.12) 

<a\ /2 \\g\\ L J\b(r>)il)\\ L2{ oy 

1/2 

We have taken into account that ||£>(D)<£>|| L2 ( ) < a 1 / HD^H^^) which can 
be checked as follows. Extend if G Hq(0; C n ) by zero to M d \ 0, keeping the 
same notation <p. Then ip G H 1 (M. d ;C n ). Using the Fourier transformation 
and the upper inequality (1.3), we obtain 

\\b(T>)<p\\l 2io) = \\b(T>)<p\\l 2{Rd) = Jjb(t)g,(t)\ 2 dt 

< «i J l£| 2 |£(0| 2 d£ = aill D( ^lli 2 (Rd) = ai||D^lli 2 (o) 
Next, by (1.2) and (2.6), 

d 

\\b(D)if)\\ L2{0) <a{ /2 J2\\Drt\\L 2{ 0) <c^ 2 d 1/2 \mm { oy (4-14) 
l=i 

From (4.12) and (4.14) it follows that 

\\A £ if)\\ H -i {0) < aidWWghJMmioy (4-15) 
Applying Lemma 4.1 to the problem (4.11), we obtain 

||r £ - if>\\m{0) < C\\A £ if)\\ H -i {0 )- (4.16) 
Now, (4.15) and (4.16) imply (4.10). • 

Remark 4.4. The statements of Lemma 4.1 and Corollary 4.2 remain true 
in any bounded domain O C M. d (without the assumption that dO G C 2 ). 
The same is true for Lemma 4.3 if the problem (4.9) is understood as the 
identity 

r (/(x)6(D)r e ,6(D)r 7 )dx = 0, Vr? G H^(0;C n ), 



L 



to 

and relation r e - if) G H%(0; C n ). 
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4.3. The "homogenized" problem. In L2(0;C n ), we consider the self- 
adjoint operator A% generated by the quadratic form 



L 



(g°b(D)u, b(D)u) dx, u G H^(0; C r 

o 



Here g° is the effective matrix defined by (1.6). Applying Corollary 4.2 with 
g e replaced by g° and taking (1.10) into account, we see that the operator 
(A ^ 1 is continuous from L 2 (0; C n ) to H^(0; C n ), and 

\\(-A°d) 1 ||l 2 (0;C™)^// 1 (C;C") < C, (4-17) 

where the constant C is defined in Lemma 4.1. Note that this fact is valid 
in any bounded domain O C R d (without the assumption that dO G C 2 ). 
Let uo G -ffo(0;C n ) be the generalized solution of the Dirichlet problem 

6(D)V6(D)u (x) = F(x), x G O; u | ao = 0, (4.18) 

where F G L 2 (0;C n ). Then u = (^) _1 F. 

Since dO G C 2 , for the solution uo of the problem (4.18) we have uo G 
H%(0;C n ) nff 2 (0;C n ), and 

ll u o ||h 2 (o ; c™) < c||F||t 2 (C;C™)- (4-19) 

Here the constant c depends only on ao, ai, H^Hloo) Ib -1 !!-^^) an d the 
domain O. To justify these properties, it suffices to note that the operator 
b(D)*g°b(D) is a strongly elliptic matrix DO and to apply the "additional 
smoothness" theorems for solutions of strongly elliptic systems (see, e. g., 
[McL, Chapter 4]). 

It follows that the operator A° D is given by the differential expression 
b(B)*g°b(D) on the domain H^(0;C n ) n H 2 (0;C n ), and that the inverse 
operator satisfies the estimate 

ll(«^£>) _1 ||L 2 (0;C")^if 2 (0;C™) < C. (4.20) 

Below we shall see that the solution u e of the problem (4.2) converges in 
L2(0; C n ) to the solution uo of the "homogenized" problem (4.18), as e — > 0. 
Our main goal is to find approximation for u e in the norm of i^ 1 (C;C ri ); 
for this, it is necessary to take the first order corrector into account. 

§5. Auxiliary statements 

In this section, we prove several auxiliary statements needed for further 
considerations. 

Lemma 5.1. Let O C M. d be a bounded domain of class C 1 . Denote 
B £ = {xeO: dist {x, dO} < e}. Then there exists a number eq G (0, 1] 
depending on the domain O such that for any u G H l (0) we have 

\ |u| 2 dx < /3e||u|| H i( )||u||£, 2 ( ), < e < e Q . (5.1) 

The constant j3 = P(0) depends only on the domain O. 
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Proof. Let us start with a model problem in the semiball T>q = 
{x G M. d : |x| < 1, Xd > 0}. For points x G M. d we write x = (x', Xd), where 
x' = (xi, . . . , Xd-i)- Introduce the following notation: 

V t = {x G R d : |x| < 1, x d > t}, S( = {x£ dV t : x d = t}, < t < e; 

T £ = {x e K d : |x| < 1, < x d < e}, S = {x G <9£> : |x| = 1}. 

Assume that u G H 1 (Vq) and u = on S. Let < i < e. Using the Green 
formula in the domain £>t, we have: 



it dx'dxd 

V t dx d 



[ \u\ 2 dx'- f 



n — — dxldxd- 

V t OXd 



Hence, 



f |n(x',t)| 2 dx' < f 2 
JSt JVt 



du 



dx d 



\u\ dx 



< 2 



du 



dx d 




Integrating over t G (0, e), we obtain 

du 



It, 



I 2 cbc < 2e 



5x d 



1/2 



dx 



(X ' ' 



: dx 



1/2 



Estimate (5.1) in the case of a bounded domain O of class C 1 is deduced 
from here in a standard way with the help of local maps, diffeomorphisms 
rectifying the boundary, and the partition of unity. Herewith, we take into 
account that the space H 1 is invariant with respect to diffeomorphisms of 
class C 1 . The number s$ must be such that the set B £o can be covered 
by a finite number of open sets admitting diffeomorphisms rectifying the 
boundary. Thus, the number eo depends only on the domain O. • 

Next statement is a direct consequence of Lemma 5.1. 
Lemma 5.2. Let O Cl' i be a bounded domain of class C 1 . Denote (dO) e = 
{x£l d : dist{x,<90} < e}. Let e x G (0, 1] be such that the set (dO) £l can 
be covered by a finite number of open sets admitting diffeomorphisms of class 
C 1 rectifying the boundary dO. Then for any u G H 1 (M. d ) we have 



I(d 



' dx < P°e\\u\\ H i 



u 



L2(ffi d )' 



< £ < £\. 



(5.2) 



The constant (3° = fi (O) depends only on the domain O. 
Proof. We apply Lemma 5.1 in the domain O and in the domain B\0, 
where B is some open ball containing O U (dO) £l . Then (5.2) is true with 
0° = max{/3(0),/3(S\0)}. • 

The following statement is similar to Lemma 2.6 from [ZhPas]. 

Lemma 5.3. Let S £ be the operator (3.1). Suppose that the domain O 
and the number e\ satisfy the assumptions of Lemma 5.2. Assume that 
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/(x) is a I '-periodic function in M. d such that f G L2(f2). Then for any 
u G ^(R^e 71 ) we kre 

|/ £ (x)| 2 |(,S £ u)(x)| 2 dx < /3*e|nr 1 ||/||i 2( n)||u|| H i (R d.cm)||u|| £a(R d. C m ) , 

e 

< £ < £ 2 , 
(5.3) 

tu/iere £2 = £i(l + n) \ /?* = + n), 2ri = diamfi. 

Proof. From (3.1), by the Cauchy inequality and the change of variables, 

we obtain 

/ |.f (x)| 2 |(S £ u)(x)| 2 dx< / dxl/^x)! 2 / |u(x-£z)| 2 dz 

<\nr' [ dy [ dz|/(£-V + z)| 2 |u(y)| 2 




u(y)| 2 dy. 



Here £ = £(1 + r\). Applying Lemma 5.2, we arrive at (5.3). • 
§6. Results in the case of bounded A 

6.1. We start with the case where Condition 1.9 is satisfied. Denote 

K° D (e) = [A'l&PX^)- 1 . (6.1) 

By (4.20), the operator b(D)(A° D )^ 1 is a continuous mapping of L2(0;C n ) 
into H 1 (0;C m ). Under Condition 1.9 the operator [A £ ] of multiplica- 
tion by the matrix- valued function A e (x) is continuous from // 1 (0;C m ) 
to H l (0;C n ). This easily follows from Corollary 2.4. Consequently, the 
operator (6.1) is continuous from L2(0;C n ) to H l {0]£- n ). 

Let u £ be the solution of the problem (4.2), and let uo be the solution of 
the problem (4.18). The "first order approximation" of u e is given by 

v e = u + £A £ 6(D)u = {A%T l ¥ + eK° D (e)F. (6.2) 
The following theorem is our main result in the case where A G L^. 
Theorem 6.1. Suppose that O C R d is a bounded domain of class C 2 . 
Let g(x) and b(D) satisfy the assumptions of Subsection 1.2. Let u e be 
the solution of the problem (4.2), and let uq be the solution of the problem 
(4.18) with F G L2(C;C n ). Suppose that A(x) is the T-periodic solution 
of the problem (1.5) and Condition 1.9 is satisfied. Let v e be the function 
defined by (6.2). Then there exists a number e\ G (0,1] depending on the 
domain O such that we have 

ll"e - V e || H i(0;C") < Co£ 1/2 ||F||i 2 ( 0;C "), < £ < £1, (6.3) 

or, in operator terms, 

Ud)s ~ (Ad)" 1 ~ eKl{e)\\ L2{0 . Cn) ^ m{0 . Cn) < Cos 1 ' 2 , < e < £l . 
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The flux p e := g £ b(D)u £ admits the following approximation 

||p e - g £ b{T>)u Q \\ L2io . Cm) < C' Q E l / 2 \\F\\ L2{0 . Cn) , < e < e u (6.4) 

where g(x) := g(x)(6(D)A(x) + l m ). T/ie constants Co, Cg depend only on 
m, d, ao, «i, IbllLoo, ||5 _1 ||Loo5 the parameters of the lattice T, the norm 
l|A|Uoo> anc ^ ^ e domain O. 

Recall that some sufficient conditions under which Condition 1.9 is sat- 
isfied are given above in Proposition 1.11. In particular, the statements of 
Theorem 6.1 are true for all operators of the form (1.1) in dimension d < 2, 
and also for the scalar elliptic operator A £ = — divg £ (x)V in arbitrary di- 
mension, where g(x) is a matrix with real entries. 

Roughening the result of Theorem 6.1, we arrive at the following corollary. 

Corollary 6.2. Under the assumptions of Theorem 6.1 we have 

l|u e - u ||l 2 (O;C") < C , o£ 1/2 ||F|| L2 ( C , ;C n), <£<£i, (6.5) 
or, in operator terms, 

\\Ad] e - {^Dy 1 \\L2{0]<C n )->L 2 {0]<C n ) < CqE 1 / 2 , < E < E\. 
~ — 1 / 2 1 1 " 

Here Co = Co + Ca 1 || , where C is defined in Lemma 4.1. 
Proof. From (6.2) and (6.3) it follows that 

K - u || L2(o) < C e 1/2 ||F|| L2(o) + e||A £ 6(D)u || L 2 (0)i < £ < £i. (6.6) 
Under Condition 1.9 we have: 

||A £ KD)u || L2( a) < ||A|| Loo ||6(D)uo||l 2(0 ). (6.7) 
Similarly to (4.13), 

||KD)uo|| L2(0) <a; /2 ||Du || L2( a). (6.8) 
Combining (6.7) and (6.8) and taking (4.17) into account, we obtain 

||A^(D)u || L2(o) < a; /2 ||A|| L J|Du || L2(o) < Ca; /2 ||A|| Loo ||F|| La(0) . 

Together with (6.6) this implies (6.5). • 

Now we distinguish the special cases. Next statement follows from The- 
orem 6.1 and Propositions 1.2 and 1.3. 

Proposition 6.3. 1°. If g° = g, i. e., relations (1.8) are satisfied, then 
A = and K^{e) = 0. In this case we have 

ll u e — "oIIh^OjC") — C , 0£ 1 ^ 2 ||F|| L2 ( C ). C n), < £ < E\. 

2°. If g° = g, i. e., relations (1.9) are satisfied, then g = g°. In this case 
we have 

||p e - 5 6(D)u || L2( o ; c-) < C^ 1/2 ||F|| L2(0 . C „), < £ < £i. 

6.2. The proof of Theorem 6.1 relies on the results for homogenization 
problem in M. d (Theorems 1.6 and 1.10) and on the tricks suggested in [Zh2, 
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ZhPas] that allow one to carry over such results to the case of a bounded 
domain. 

Let us fix a linear continuous extension operator 

Po : H 2 (0;C n ) -> H 2 (R d ;C n ), (6.9) 
and put uo = Po u o- Then 

l|uo||// 2 (R d ;C' 1 ) ^ CoH u o||h 2 (C';C™)> (6.10) 

where Co is the norm of the operator (6.9). Denote 

v«(x) = uo(x) + eA £ (x)6(D)S (x). (6.11) 

Then v e = v^ 1 ' \q. 

The following statement is proved with the help of Theorems 1.6 and 1.10. 

Lemma 6.4. Let uq be the solution of the problem (4.18), and let v £ be the 
function defined by (6.2). Then for < e < 1 we have 

\\A £ V £ - yl Uo||/f-i(o ; C") < C4£||uo||if2( 0;C n). (6-12) 

The constant C4 depends only on m, d, Qo, «i, H^Hl,*,, lls" 1 !!-^' ^ e Pa- 
rameters of the lattice T, the norm HAH^, and the domain O. 

Proof. The required estimate in the case of a bounded domain is deduced 
from the similar inequality in W 1 . Let Vg 1 "* be defined by (6.11). We check 
that 

\\A £ v^ - _4°5 || H -i (R d. C n) < C^e\\uo\\ H 2 {Rd . cn) , < e < 1. (6.13) 
Clearly, 

F := ^°S + S € L 2 (R d ; C n ). (6.14) 
By the Fourier transformation and (1.3), (1.10), we obtain 

\\nl m = [ |(6(OV6(o + i)fio«)l 2 de 

< [ (aii5 ||«e| 2 + 1) 2 |Q (€)| 2 ^ < (max{ Ql || 5 || ioo ,l}) 2 ||S ||^ 2(IRd) . 

(6-15) 

Here uo(£) is the Fourier-image of the function Uo(x). 

Let s £ € i2' 1 (M <i ; C n ) be the generalized solution of the equation 

As £ + s e = F. (6.16) 

From Theorems 1.6 and 1.10 it follows that 

ll s e - uo||L 2 (Rd ; c™) < CiellFllLa^C"). < e < 1, (6.17) 

ll S £ - V^HiJl^d.Cn) < C3£||F|| L2(M d ; cn), 0<e<l. (6.18) 
By (6.14) and (6.16), 

A £ w^ - ^°S = A £ (vW - s £ ) + A £ s £ - A°u = A £ (v £ ^ - s £ ) - (s £ - S ). 



2d 
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Hence, 

\\A £ vP -A°u \\ H -i {Rd) < \\A £ (vP-s £ )\\ H -i (Rd) + \\s £ -u \\ H -i {Rd) . (6.19) 
Next, taking (1.3) into account, we obtain 



(/6(D)(v| 1) -s £ ),6(D)r 7 



Li 



< "iIMIlJIvW - s £ \\ H i 

Combining this with (6.17)-(6.19), we see that 

\\A £ vP-A u \\ H - 1{Rd) < (Ci+Csaill^ll^JellFH^^), < e < 1. (6.20) 

Now, (6.15) and (6.20) imply (6.13) with the constant 

C A = (Ci +C 3 ai||5(|| Loo )max{ai|| fi r|| Loo ,l}. 

Returning to the case of a bounded domain, note that if f £ -ff _1 (0;C n ) 
and f € H' 1 ^; C n ) are such that f \ = f, then 

\f n (f,<p)dx\ Ld(f,<P)dx 

ll f llH-i(o) = su p ii II = su p — 

o^(p&cg°(o) W^PWmio) o^ v ecg°(o) \\V > \\H 1 («. d ) 

f Rd (?,<p)dx 
— SU P — ii — ii 

0^ipeCg°(R d ) \\ ( P\\H 1 (R d ) 

Hence, 

\\As% - A uo\\ H -i( ) < || Av^ - A°u \\ H -i( Rd y 
Together with (6.13) and (6.10) this yields 

||Av e — «4°Uo||if-i(0) < C4£||u ||^2( R d) < C±Co£\\xi \\ H 2( O y 
Thus, inequality (6.12) holds with C4 = C\Cq- • 

6.3. The fist order approximation v £ of the solution u £ defined by (6.2) does 
not satisfy the Dirichlet condition on dO. We consider the "discrepancy" 
w £ which is the generalized solution of the problem 

A £ w £ = in O, w £ \ao = v £ \ do = eA e 6(D)u |aci. (6.21) 

Here the equation is understood in the weak sense: the function w £ € 
H x {0\C n ) satisfies the identity 

/ (/(x)6(D)w £ ,6(D)77) dx = 0, Vrj € Hq(0; C n ). 
Jo 

The boundary condition in (6.21) is understood in the sense of the trace theo- 
rem: under the assumptions of Theorem 6.1 one has A £ 6(D)uo 6 H 1 (0;C™), 
whence A £ 6(D)u | 9o G # 1/2 (d0; C n ). 



IflljT-i 
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By (4.2) and (4.18), .4 £ (u £ -v £ ) = _4°u -.4 £ v £ . Consequently, by (6.21), 
the function u £ — v £ + w £ is the solution of the following Dirichlet problem 

A £ (u £ -v £ + w £ ) = A°u - A £ v £ in O, (u £ - v £ + w £ )|ae> = 0. 

The right-hand side in the equation belongs to H~ l (0; C n ). Then, applying 
Lemmas 4.1 and 6.4, for < e < 1 we obtain 

||u £ - V £ + W £ \\ H i^ ;C n ) < CH^Uo - «4eV e ||fi--:i(0;C») < CC4£\\u \\ H 2( O . C ny 

Together with (4.19) this implies that 

||u e - v £ + w e || H i (0;Cn ) < CC±c e||F|| i2 ( 0;C n), < e < 1. (6.22) 

Therefore, the proof of estimate (6.3) from Theorem 6.1 is reduced to esti- 
mating of w £ in H l (0; C n ). 

Assume that < e < £i, where the number E\ € (0, 1] is defined in Lemma 
5.2. Fix a smooth cut-off function £ (x) in M. d supported in the e-vicinity 
of the boundary dO and such that 

6 £ € C?(R d ), supple (dO) e , O<0 £ (x)<l, 
^£( x )|ao = 1) e | V6' e (x) | < re = const. 

Consider the following function in 

e (x) = e# £ (x)A £ (x)6(D)5 (x). (6.24) 

Then £ G H l (R d ;C n ) and e | 9O = eA £ &(D)u |ae>. The problem (6.21) can 
be rewritten as: -4 £ w £ = in O, w £ \do = </>e\dO- Applying Lemma 4.3, we 
obtain 

l|w £ ||iji(o ; C") < lo\\^ E \\m(0;C n )- (6.25) 

Thus, the proof of the required estimate for the norm of w £ in ff 1 (0; C n ) is 
reduced to the next statement. 

Lemma 6.5. Suppose that the assumptions of Theorem 6.1 are satisfied. 
Assume that < e < e±, where the number e± € (0, 1] is defined in Lemma 
5.2. Let 4> £ be the function defined in accordance with (6.23), (6.24). Then 
we have 

He\\m(0;C") < C 5 e 1/2 \\F\\ L2(0 . cn) , < e < £l . (6.26) 

The constant C5 depends only on m, d, ao, a±, H^Hl^, lis -1 Hl^ , the norm 
ll-^-IUoo) an d the domain O. 

Proof. The norm of <f> £ in L2(0; C n ) is estimated with the help of Condition 
1.9 and relations (4.17), (6.8), and (6.23): 

H £ \\ L2{0 ) <e||A £ 6(D)u || L2(o) <e||A|| Loo ||6(D)u || L2 ( O ) 
< eai /2 ||A|| Loo ||uo|| H i(o) < eCai /2 ||A|| Loo ||F|| £2(0) . 
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Consider the derivatives 

d~ct> £ de £ 



dxj 



£ |^A £ &(D)G + 9 £ ( J^-) £ fe(D)S + ee £ A%b(B)d 3 u ), 

j = l,...,d. 

Then 

||D0 £ ||i 2(o) < 3e 2 J |W £ | 2 |A £ 6(D)u | 2 dx + 3 J l(DA) £ |> e 6(D)5 | 2 

d ,. 

+ 3e 2 V / \O e \ 2 \\ e b(p)Djvm\ 2 dx.. 
Pi Jo 

(6.28) 

Denote the terms in the right-hand side of (6.28) by X l5 1 2 , and Z 3 , respec- 
tively. 

It is easy to estimate X 3 . By (6.23), Condition 1.9, and (1.2), (2.6), we 
obtain 

d 

||^A £ 6(D)^u || 2 2(o) < IIAH^aid^llA^uolli,^). 

i=i 

Together with (4.19) this yields 

Z 3 < Se^lAH^aidHuoll^^) < 73£ 2 ||F||| 2(0) , (6.29) 

where 73 = 3c 2 o;id||A||| . 

In order to estimate the first term in the right-hand side of (6.28), we 
apply (6.23), Condition 1.9 and Lemma 5.1. We have 

Zi < S^IIAHl^ / |6(D)u | 2 dx 

J B s 

<3 K 2 ||A|| 2 oo /3 e ||6(D)uo||H 1 (o)l|fe(D)u || L2( o ) . 

Using (4.17), (4.19), (6.8), and the estimate 

d 

||6(D)uo||hi(o) = E 6 iAuo|| Td 1/2 \\MH> { 0), 
1=1 

we arrive at the inequality 

X 1 < fe^HAII^^aid^lluoUHHoKll^co) < 7ie||F||£ 2(0) , (6.30) 
where 71 = 3cCk 2 ||A|| 2 paid 1 ' 2 . 



— " e 6(D)uo) 



dx 



2 



dx. 
(6.31) 
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It remains to consider the second term in the right-hand side of (6.28). 
By Corollary 2.4, 

Z 2 <3 f \(DA) £ \ 2 \9 e b(D)u \ 2 dx 

~ d 

< 3A ||^6(D)5 ||| 2(]Rd) + S^IIAHi^^ 2 / 
Since 

r) f)f) B 

_ ( « D)ao) = _£ KD)ao + e ,_ (KD)ao) , 
then, by (6.23) and (6.31), we have 

X 2 < 3 (ft + 2/3 2 ||A||! oo k 2 ) / |6(D)S | 2 dx 

J(dO) € 

,. d 

+ 6/3 2 ||A||| oo e 2 / V|6(D)Au | 2 dx. 

Combining this with Lemma 5.2 and condition (1.3), we obtain 
X 2 < 3 (ft + 2/3 2 ||A||! oo k 2 ) /3 e ||6(D)5o||^i (Md) ||fe(D)S |lL 2 
+ 6ft||A||l oo e 2 ai||5o|| 2 ^ 2(Rd) 

< 3e (ft + 2/3 2 ||A||| oo k 2 ) /^Vlluoll^QR^IIuolltfi^d) 

+ 6e 2 /3 2 ||A||| oo ai||5 || 2 ^2 (Rd) . 

Taking (4.19) and (6.10) into account, from (6.32) we deduce that 

X 2 < 7 2£||F||| 2(0) , (6.33) 

where 72 = 3(cC ) 2 ((ft + 2ft||A|| 2 oo K 2 ) /?V + 2ftai||A|| 2 J. 
Now, relations (6.28)-(6.30) and (6.33) imply that 

\\n^ e \\l 2(0) <l 1 +l 2 +l 3 <e( 7 i + 72 + 73)||F|| 2 2{0) , <e< £l . (6.34) 

Finally, (6.27) and (6.34) yield estimate (6.26) with 

1/2 



(6.32) 



C 5 = (6 2 ai \\A\\ 2 Loo + 7i + 72 + 7s) ' 



Now it is easy to complete the proof of Theorem 6.1. From (6.22), 
(6.25), and (6.26) it follows that 

Hue - VeIIh^OjC") - ^^4ce||F|| Z/2 ( C i ;C n)-|-7oC'5e 1 / 2 ||F|| L2 ( C ). C n), < e < e±, 

which implies (6.3) with Co = CC4C + 70C5. 

It remains to check (6.4). From (6.3), (1.2), and (2.6) it follows that 

|| Pe -5 £ 6(D)u -/6(D)(eA £ 6(D)u )|| L2(o) 
< \\9\\L^ 2 d l l 2 C Q e l l 2 \\¥\\ L2{0) , 0<e<e 1 . 
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From (1.2) and the definition of the matrix g it is seen that 

d 

/6(D)u +/6(D)(eA £ 6(D)u ) = 5 £ &(D)u +e 5 £ ^6 i A £ 6(D)Auo. (6.36) 

l=i 

Applying Condition 1.9 and relations (1.2), (2.6), and (4.19), we obtain 

d 

\\eg e ^A £ 6(D)Au || L2(o) < e\\g\\ Lgo ||A|| £oo aidc||F|| L2(0) . (6.37) 
i=i 

Now, relations (6.35)-(6.37) imply (6.4) with the constant Cg = 
h\\L oa oi] /2 d 1 l 2 Co + |b||Lj|A|| Loo aic?c. • 

§7. Results in the general case 

7.1. Now we refuse the assumption that A(x) is bounded. Then we need to 
include a smoothing operator in the corrector. 

Let Pq be the extension operator (6.9), and let S £ be the operator smooth- 
ing in Steklov's sense defined by (3.1). By Rq we denote the operator of 
restriction of functions in W 1 to the domain O. We put 

K D (e) = i? [A £ ]S £ 6(D)P (^)- 1 . (7.1) 

The operator b(D)Po(A^ ) )~ 1 is a continuous mapping of L2(0;C™) into 
H 1 (M. d ;C m ). As has been mentioned in Subsection 3.2, the operator [A e ]S £ 
is continuous from H 1 (M d ;C m ) to if 1 (IR d ; C"). Consequently, the operator 
(7.1) is continuous from L 2 (0;C") to H\0;C n ). 

Let u £ be the solution of the problem (4.2), and let uo be the solution of 
the problem (4.18). As above, we denote So = Po u o- We put 

vf(x) = uo(x) + £ A £ (x)(5 e 6(D)5 )(x), 

and v £ := v £ 2 ^|e>. Then 

v £ = (A° D )- 1 F + eK D (e)F. (7.2) 
The following theorem is our main result in the general case. 

Theorem 7.1. Suppose that O C M. d is a bounded domain of class C 2 . 
Let g(x) and 6(D) satisfy the assumptions of Subsection 1.2. Let u e be 
the solution of the problem (4.2), and let u$ be the solution of the problem 
(4.18) with F G L 2 (0;C n ). Let v £ be the function defined by (7.1), (7.2). 
Then there exists a number e 2 € (0, 1] depending on the domain O and the 
lattice T such that we have 

K - v e || H i (0 . c ™) < Ce 1/2 ||F|| L2(0;Cn) , 0<e<£ 2 , (7.3) 

or, in operator terms, 

W^D,e ~ (^Dr 1 - £K d(£)\\l 2 (0;C")^HI{0;C") < Ce 1/2 . 

The flux p £ := g £ b(D)u £ admits the following approximation 

||p £ - g £ S £ b{T))^\\ L2{0 . €m) < C'e^llFlli^c-cn), < e < e 2 , (7.4) 
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where g(x) := g(x)(6(D)A(x) + l m ). The constants C, C depend only on m, 
d, ao, cxi, IblUoo, ||5 _1 ||Loo5 the parameters of the lattice T, and the domain 
O. 

Roughening the result of Theorem 7.1, we arrive at the following corollary. 

Corollary 7.2. Under the assumptions of Theorem 7.1 for < e < £2 we 

have 

\\ u e ~ Uo||l 2 (0;C") < Ce 1/2 ||F ||l 2 (0;C") 1 ( 7 - 5 ) 
or, in operator terms, 

W^D,e ~ (-^d)~ 1 |Il2(C;C")^l 2 (C;C") < Ct 1 / 2 . 

T/ie constant C is given by 

n r* 1 n ^ l/2/o \-l - 1 / 2 1/2 M 1,1/2,, -i M l/2 
C = C + C cm' (2r ) a a/ HsH/JIsf || Loo , 

where c is the constant from (4.19), Co is the norm of the extension operator 
Po, and ro is the radius of the ball inscribed in closfi. 
Proof. From (7.2) and (7.3) it follows that 

K " uo||l 2( o) < Ce 1/2 ||F|| L2(0 ) + e\\A £ S e b(D)u \\ 
By (3.9) and (1.3), 

\\A £ S £ b(B)u \\ L2{o) < \\A s S £ b(B)u \\ L2{Rd) 

1/2 ~ ( 7 - 7 ) 
< Af||6(D)uo||£ 2 (R<i) < Ma{ ||uo|| H i(Rd). 

Taking (4.19) and (6.10) into account, we obtain 

l|2o||tfi(R d ) - l|2o||H 2 (R«i) < Co||uo||// 2 (0) < Coc||F II L 2 (0)- (7.8) 

Now, from (7.6)-(7.8) it follows that 

K " u oIIl 2 (o) < C£ 1/2 ||F|| L2(0) + eMa\ /2 C c\\nL 2 (0)- 
Recalling the expression for M (see (3.8)), we arrive at (7.5). • 

7.2. Let us start the proof of Theorem 7.1. The following statement is 
similar to Lemma 6.4. 

Lemma 7.3. Let uq be the solution of the problem (4.18), and let v £ be the 
function defined by (7.1), (7.2). Then for < e < 1 we have 

u o||H- 1 (C';C n ) — C6 £ || u o||H 2 (£>;C n )- 
Here the constant C% is given by 

C 6 = Ce>(Ci + (7 2 ai||c/||L 00 )max{ai||5r||L 00 ,l} 

and depends only on m, d, ao, ct±, H^Hl^, Wd" 1 1 the parameters of the 
lattice r, and the domain O. 

Proof. Lemma 7.3 can be proved by analogy with the proof of Lemma 
6.4. The only difference is that one should apply Theorem 3.3 instead of 
Theorem 1.10. • 
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Next, by analogy with the proof of Theorem 6.1, we consider the "dis- 
crepancy" w £ G H l (0; C n ) which is the generalized solution of the problem 

,4 e w e = mO, w £ \ do = v £ \ do = eA £ (S £ b(B)u )\ do . (7.9) 

The equation in (7.9) is understood in the weak sense, and the boundary 
condition in the sense of the trace theorem. It should be taken into account 
that A £ (S £ b(B)u ) e H l (Q-C n ). 

By (4.2), (4.18), and (7.9), the function u e — v e + w e is the solution of 
the following problem 

i e (u £ - v £ + w E ) = i°u - i £ v £ in O, (u e - v e + w e )\ do = 0. 

Applying Lemmas 4.1 and 7.3, for < e < 1 we obtain 

ll u £ — v e + w e\\m(0;C n ) < C7||^4°Uo — -^e v e\\H- 1 {0;C n ) < CCQ£\\u \\ H 2^ o . €n y 

Together with (4.19) this implies that 

||u e - v e + w £ \\ H i {0;Cn) < CC 6 ce||F|| L2(c , ;C u), < e < 1. (7.10) 

7.3. By (7.10), the proof of estimate (7.3) is reduced to estimating of the H 1 - 
norm of w £ . As in Subsection 6.3, we fix a cut-off function 6> e (x) satisfying 
conditions (6.23). We assume that < e < £2, where the number £2 £ (0, 1] 
is defined in Lemma 5.3. Consider the following function in R. d : 

e (x) = £0 £ (x)A £ (x)(S £ 6(D)uo)(x). (7.11) 

Similarly to (6.25), by Lemma 4.3, we have 

ll w £||// 1 (C;C n ) - 7o||0 e ||#i(O;C n )- ( 7 -12) 

Thus, the problem is reduced to the proof of the following statement. 

Lemma 7.4. Suppose that the assumptions of Theorem 7.1 are satisfied. 
Let < e < £2, where the number £2 G (0, 1] is defined in Lemma 5.3. Let 
4> £ be the function defined in accordance with (6.23), (7.11). Then we have 

H^ellif^OjC") - C 7£ 1/2 \\F\\ L2 ( ;C n )' < £ < £ 2 . (7.13) 

The constant Cj depends only on m, d, ao, a±, H^Hl^,, ||<7 -1 ||l 00) the pa- 
rameters of the lattice T, and the domain O. 

Proof. We start with the estimate for the norm of the function (7.11) in 
L 2 (C;C n ). From (1.3), (3.9), (6.23), and (7.8) it follows that 

\\4> £ \\l 2 (o) < 4 A£ (S £ b(D)u )\\ L2{Rd) < £Ma} /2 ||So|| H i (Rd) 



< £Ma{ /2 C c\\F\\ L2i0) . 
Consider the derivatives 

i>0 - - £^A £ (S £ b(B)u ) + 9 £ ( (S £ b(D)u Q ) 



OXj OXj \ OXj J 

+ £^A £ (5 £ 6(D)^u ), j = l,...,d. 
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Then 

HD^IlL(o) <3e 2 ^|V^| 2 |A £ (S' £ 6(D)So)| 2 dx 

+ 3/ \(T>A) £ \ 2 \e £ (S £ b(-D)u )\ 2 dx + 3e 2 Y [ % | 2 |A £ (5 e 6(D)DjUo) | 2 dx. 

(7.15) 

The summands in the right-hand side of (7.15) are denoted by J\, J2, and 
Jz-, respectively. 

It is easy to estimate J3. From (1.3), (3.9), and (6.23) it follows that 

d 

Js < 3e 2 l|A £ (5 , £ &(D) J D,5 )|| 2 2(Rd) < 3e 2 M 2 a 1 ||S ||^ (Rd) . 
3=1 

Combining this with (4.19) and (6.10), we obtain 

Jz<%e 2 \\nl 2i oy (7-16) 

where 73 = 3M 2 ai(C c) 2 . 

The fist term in the right-hand side of (7.15) is estimated with the help 
of (6.23) and Lemma 5.3. For < e < £2 we have 



Ji < 3k 2 f \A £ (S £ b(D)u ) | 2 dx 



<{dO) e 

< 3/c 2 ^*e|n| _1 ||A||| 2( ^ ) ||6(ID)uo||i?i(ffid ) i|&(D)uo|| jL2(]R rf ) . 

Combining this with (1.3), (4.19), (6.10), and estimate (2.14), we arrive at 
the inequality 

(7.17) 

where 71 = 3k 2 ^(C o c) 2 m{2r )- 2 a Q 1 a 1 \\g\\ Loo \\g~ l \\ Lgo . 

It remains to consider the second term in the right-hand side of (7.15). 
By (6.23), 

J 2 <sf |(DA) £ | 2 |S £ 6(D)u | 2 dx. 

J(dO) £ 

By Lemma 5.3, for < e < £2 we have 

J2<3/3^|O|- 1 ||DA|| 2 2(n) ||fe(D)S || 
Together with (1.3), (2.15), (4.19), and (6.10) this implies that 

J2 <72£||F|| 2 2(0) , 0<£<£ 2 , (7.18) 

where 72 = 3/3*(C C ic) 2 ma " 1 ai||5(||L oo Hs^llw 
Finally, relations (7.15)-(7.18) yield 

||D0 e || 2 2(o) <Ji + J 2 + Jz< (71 + 72 + 73)£||F|| 2 2(0) , < £ < £ 2 . 
Combining this with (7.14), we obtain (7.13) with 

C 7 = {M 2 ai {C c) 2 + 71 + 72 + 7s) 1/2 - • 
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Now, it is easy to complete the proof of Theorem 7.1. From (7.10), 
(7.12), and (7.13) it follows that 



\ u e-^e\\m(0;C n ) < CC e ce\\F\\ L2 ^ . cn ) +7 CVe 1 / 2 ||F|| L2 ( 0;C n) , < e < 



£2- 



This implies (7.3) with C = CC 6 c + j C 7 . 

It remains to check (7.4). Taking (1.2) and (2.6) into account, from (7.3) 
we obtain 

||p e - <? £ 6(D)u - g £ b(D)(sA £ S £ b(D)u )\\ L2(o) 

(7 19) 

^IblUco^V/W^HFiu^o), o< £ < £2 . 

By Proposition 3.1 and relations (1.3), (4.19), and (6.10), we conclude that 

||<7 £ 6(D)u - 5 e 5 e 6(D)5o|| L2(0 ) < \\g\\ Loo \\b(D)u - S £ b(B)u \\ L2{Rd) 

< e r i||5llLocai /2 ||5olb2 (Md) < eri||5|| Loo ai /2 C cl|F|| L2(0) . 

(7.20) 

From (1.2) and the definition of the matrix g it is seen that 
g E S £ b(D)u + g £ b(T>)(eA £ S £ b(T>)u ) 

= g £ S £ b(B)u + eg £ £ &,A £ ^(D)A5o. 
i=i 

Taking (1.3), (2.6), (3.9), (4.19), and (6.10) into account, we obtain 
d 

H 2 (R d ) / voo \ 

1=1 ~ (7.22) 

<e\\g\\ Loc M ai d l l 2 Coc\\nL 2 (o)- 

Now, relations (7.19)-(7.22) imply (7.4) with the constant C = 
Uk^d^C + hW^Cocinal 12 + Ma x d^). . 
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